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Abstract 



The new theoretical input to the analysis of the experimental data of the CCFR col- 
laboration for F$ structure function of uN deep inelastic scattering is considered. This 
input comes from the next-to-next-to-leading order corrections to the anomalous dimen- 
sions of the Mellin moments of the F% structure function and N 3 LO corrections to the 
related coefficient funtions. The QCD scale parameter A^^- is extracted from higher-twist 
independent fits. The results obtained demonstrate the minimization of the influence of 
perturbative QCD contributions to the value of A^-. 



1 Introduction 



One of the most important current problems of symbolic perturbative QCD studies is the 
analytical evaluation of the next-to-next-to-leading order (NNLO) QCD corrections to the 
kernels of the DGLAP equations for different structure functions of the deep-inelastic 
scattering (DIS) process. In this note we will apply the related information for the fixation of 
definite uncertainties of the NNLO analysis p], || of experimental data for F 3 structure function 
(SF) data of uN DIS, provided by the CCFR collaboration |4|] at the Fermilab Tevatron and 
present preliminary results of our improved fits which will be described elsewhere ||. 



2 Methods of analysis of DIS data 



There are several methods of analysis of the experimental data of DIS in the high orders of 
perturbation theory. The traditional method is based on the solution of the DGLAP equation, 
which in the case of the F 3 SF has the following form: 



y 



V F3 (y,A s )+f3(A & 



d\nC Fa (y,A s 
dA» 



x 



Fs[-,Q 2 

y 



(i) 



l 



where A s = a s / (47r), fidA s /dfi = (3(A S ) is the QCD /^-function and Cp 3 (y, A s ) is the coefficient 
function, defined as 

C F3 (y,A s ) = ^C F ay)(^y (2) 

and Vp 3 (z) is the DGLAP kernel, related to a non-singlet (NS) F 3 SF. The solution of Eq.(l) 
is describing the predicted by perturbative QCD violation of scaling || or automedeling 
behaviour of the DIS SFs by the logarithmically decreasing order a s -corrections. 

The coefficient function we are interested in has been known at the NNLO for quite a 
long period. The term Cj? 3 ,2(?/) was analytically calculated in Ref. ||. The results of these 
calculations were confirmed recently using a different technique. 

The kernel Vp 3 {z,a s ) is analytically known only at the NLO. However, since there exists a 
method of symbolic evaluation of multiloop corrections to the renormalization group functions 
in the MS-scheme fl(| and its realization at the FORM system, it became possible to calculate 
analytically the NNLO corrections to the n — 2,4, 6, 8, 10 Mellin moments of the NS kernel of 
the F% SF fH|| . They have the following expansion: 



f z n ~ l V NS , F2 {z,a s )dz = E7SU(^)(^) ,+1 ( 3 ) 



and are related to the anomalous dimension of NS renormalization group (RG) constants of F 2 
SFR : 

d\nZ NS ^ , ^ 
V 9^ =Tns,fMs) ■ (4) 

These results were used in the process of the fits of Refs. @, |3| of the CCFR data for the F 3 
SF with the help of the Jacobi polynomial method |13| . It allows the reconstruction of the SF 
F3 from the finite number of Mellin moments Mj i p 3 (Q 2 ) of the xF% SF: 

ax ?l 

F^ ma *{x,Q 2 ) = wJ2 e^(x)^cf(a,/3)M™J 3 (Q 2 ) (5) 

n=0 j=0 

where w = w(a,[3) = x a_1 (l — x)@, 0" ,/3 are the orthogonal Jacobi polynomials and c^ l \a,(3) 
is the combination of Euler T-functions, which is factorially increasing with increasing of N max 
and thus n. 

The expressions for Mj^^iQ 2 ) include the information about Mellin moments of the co- 
efficient function 

C n , F3 (Q 2 ) = C x n - 1 C F3 {x,a s )dx = Y J C {i) {n)(^-\ (6) 
Jo \AttJ 

where C^(n) = 1. The target mass corrections, proportional to (M^/Q 2 )M J+4j p 3 (Q 2 ), are also 
included into the fits. Therefore, the number of the Jacobi polynomials N max = 6 corresponds 

1 The method of renormalization group was originally developed in fll2|. 
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to taking into account the information about RG evolution of 10 moments, and N max = 9 
presumes that the evolution of n = 13 number of Mellin moments is considered. 

The procedure of reconstruction of F%(x, Q 2 ) from the finite number of Mellin moments and 
the related fits of the experimental data were implemented in the form of FORTRAN programs. 
The details of the fits of the CCFR data, based on RG evolution of 10 moments, are desribed 
in Refs.p, |3| (for the brief review see Ref.||14||). In the process of these analyses the following 
approximations were made: a) it was assumed that for a large enough number of moments, 
jNSF 3 ( a s) ~ 7vs , F 2 ( a s)i b) si nce the odd NNLO terms of 7yvsF 2 are explicitly unknown, they 



were fixed using the smooth interpolation procedure proposed in Ref.[|15|. It was known that 
the additional contributions, proportional to the d abc d abc structure of the colour gauge group 
SU(N C ) are starting to contribute to the coefficients of 7^s,F 3 ( a s) from the NNLO 0. In the 
process of the analysis of Refs. 0, |3[] it was assumed that they were not dominating and therefore 
were not taken into account. 



3 New inputs for the fits 

After recent explicit analytical evaluation of the NNLO coefficients of 7^5 F 3 ( a s) & t n =3,5,7 
,9,11,13 (see Ref.[|l6|) it became possible to fix this uncertainty (it is worth noting that the 
NNLO contribution to F 2 ( a s) f° r n = 12 was analytically evaluated in Ref . |TB|| also). To 
estimate the NNLO terms of 7^5 Fs (a s ) at n =4,6,8,10,12 we applied the smooth interpolation 
procedure, identical to the one used to estimate the odd NNLO terms of ^l F2 (a s ), while the 

numerical value of 7^5_p 3 (2) was fixed with the help of an extrapolation procedure, where we 
have not used the value at n — 1. The justification and more details of this procedure will be 
given elsewhere ||. 

The used numerical results of the NNLO contributions 7^5. F 3 ( n ) with and without d abc d abc - 
factors are presented in Table 1, where we marked in parenthesis the estimated even terms. 
The expressions for the NNLO contributions to the NS anomalous dimensions terms ^ NS F2 (n) 
are also given for comparison. They include the numerical results of the explicit analytical cal- 
culations of Refs. [|Tl] , [TBI , normalized to / = 4 numbers of active flavours, and the results of the 
smooth interpolation procedure, in parenthesis, applied for estimating explicitly uncalculated 
odd terms. The satisfactory agreement between the numbers in the second and third columns 
supports the assumptions a) and b) mentioned above. 

In Table 2 the numerical expressions for the coefficients of Eq.(6) for / = 4 numbers of 
active flavours are given. They include the results of explicit calculations of N 3 LO corrections 
of odd moments supplemented with the information about the coefficients of the Gross- 
Llewellyn Smith sum rule [17], defined by the n = 1 Mellin moment of the xF^ SF. The 
numbers in parenthesis are the results of the interpolation procedure. In the last column we 
present the values of C^(n), obtained with the help of the [1/1] Pade estimates approach. One 
can see that the agreement of Pade estimates with the used N 3 LO results is good in the case 



3 



n 


i NS,Fs\ J 


d abc d abc neglected in 


'vlro n ("ft) 
lNS,F3\'"/ 


'vlro n (?l) 
lNS,F2\ > 


2 


(631) 


(585) 




612.06 


3 


861.65 


836.34 




(838 93) 

\ v7 - 7 ^7 . t_7 -.7 / 


4 


H 015 37) 


fl 001 42) 




1005 82 


5 


1140.90 


1132.73 




(1135.28) 


6 


(1247) 


(1241.21) 




1242.01 


7 


1338.27 


1334.32 




(1334.65) 


8 


(1420) 


(1416.73) 




1417.45 


9 


1493.47 


1491.13 




(1492.02) 


10 


(1561) 


(1558.85) 




1559.01 


11 


1622.28 


1620.73 




(1619.83) 


12 


(1679.81) 


(1677.70) 




1678.40 



Table 1: The numerical expressions of the NNLO coefficients of anomalous dimensions of the 
n-th NS moments of the F 3 and F 2 SFs at / = 4. The numbers in parenthesis are the estimated 
results. 
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-1.78 
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4 


4.87 
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13 
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15923.91 


17353 



Table 2: The numerical expressions for the coefficients of the coefficient functions for n-th 
Mellin moments of the F 3 SF up to N 3 LO and their [1/1] Pade estimates. 
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Table 3: The results of the fits of the CCFR data for xF% SF, taking into account the NNLO 
approximation for jp^NS' The initial scale of RG evolution is Ql=20 GeV 2 . 



of the Gross-Llewellyn Smith sum rule (this fact was already known from the considerations 
of Ref.||19||). In the case of n = 2 and n > 6 moments the results are also in satisfactory 
agreement. Indeed, one should keep in mind that the difference between the results of column 
3 and 4 of Table 2 should be devided by the factor (1/4) 3 , which comes from our definition of 
expansion parameter A s = a s /(Air). Note, that starting from n > 6 the results of application 
of [0/2] Pade approximants, which in accordance with analysis of Ref. |^0] are reducing scale- 
dependence uncertainties, are even closer to the the results of the interploation procedure (for 
the comparison of the estimates, given by [1/1] and [0/2] Pade approximants in the case of 
moments of xF 3 SF see Ref. 0, while in Ref. |^T| the similar topic was analysed within the 
quantum mechanic model). For n =3,4 the interpolation method gives completely different 
results. The failure of the application of the Pade estimates approach in these cases might be 
related to the irregular sign structure of the perturbative series under consideration. 



4 Some results of the fits 



In Table 3 we present the comparison of the results of the determination of the parameter, 
made in Ref. ||, with the new ones, obtained by taking into account more definite theoretical 
information. Since NNLO corrections to the anomalous dimensions and N 3 LO contributions to 
the coefficient functions of odd moments of the xF 3 SF are now known up to n =13, it became 
possible to study the dependence of the results of the fits from the value of N max , which we 
can now vary from N max = 6 to N max = 9. It should be mentioned that for N max = 6 the new 
NNLO result and its Q$ dependence are in agreement with the results of Ref. 0. However, the 
incorporation of higher number of moments, and thus the increase of N max , make the NNLO 
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(and approximate N 3 LO ) results almost independent from the variation of Qq in the interval 
5 GeV 2 -100 GeV 2 . This is the welcome feature of including into the fits the results of the new 
analytical calculations of the NNLO corrections to anomalous dimensions and N 3 LO corrections 
to the coefficient functions of odd moments of the xF% SF [jT6| . Comparing now the central 
values of the results of the stable NLO fits of Ref.|| with the new NNLO and N 3 LO results, 
we observe the decrease of the theoretical uncertainties and, probably, the saturation of the 
predictive power of the corresponding perturbative series at the 4-loop level. More detailed 
results of our fits, including extraction of a s (M z ), its scale dependence and the information 
about the behaviour of twist-4 corrections at the NNLO and N 3 LO, in the case of N max = 9, 
will be described elsewhere 
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